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Сформулированы основные недостатки метода дифференциально-тейлоровских преобразований и осуществлен переход от степенного базиса 
к базисам ортогональных многочленов. Показано, что при осуществлении перехода к разложению по многочленам Чебышева первого рода и 
смещенным многочленам Чебышева существенно увеличивается скорость сходимости ряда. Сформулирован более универсальный алгоритм 
вычисления дискрет дифференциального спектра. Показано, что в чебышевских базисах величина дискрет спектра постоянно уменьшается 
с увеличением их номера. В этом случае можно остановить вычисление дискрет по достижению их величины требуемого малого значения, 
что нельзя сделать в степенном базисе. На численных примерах показаны преимущества перехода к разложению по многочленам Чебышева 
первого рода и смещенным многочленам Чебышева. 


1. Дифференциально-тейлоровские преобразования функций 

В тех задачах математического моделирования, в 
которых решение уравнений на заданном интервале 
изменения независимого переменного эффективно 
может быть представлено степенными рядами Тейлора, 
нашли применение дифференциально-тейлоровские 
преобразования [1]: 
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где слева от символа л г стоит прямое 
ванне оригинала х(1) в изображение Х(к) - дискреты 
дифференциального спектра, а справа - обратное 
преобразование Х(к) в х(1), к - дискретный аргумент і 
- время, Н, т - постоянные. 

Для восстановления х(1) по результатам измерений 
хДД, х(1 2 ),. . .,х(Ц) необходимо решить систему линейных 
алгебраических уравнений вида 
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из которой определяются неизвестные дискреты Х(1 ^ ^ 

Главным недостатком дифференциальных пре- 
образований является то, что при исследовании вре- 
менных характеристик реальных систем, в процессе 
восстановления оригинала по формуле (1) возникает 
необходимость возведения переменной 1; в большие 
степени. Для 1; порядка микросекунд или наносекунд 
становится проблематичным учет более десяти-пятнад- 
цати членов ряда, так как система (2) становится плохо- 
обусловленной, и ошибка вычислений при операциях с 
числами большой или малой размерности становится 
недопустимой. Попытки определить оптимальное коли- 
чество вычисленных дискрет с точки зрения получения 
наибольшей точности решения в степенном базисе 
оказались безуспешными. 

Соседние дискреты, как показали многочисленные 


расчеты, могут отличаться на несколько порядков. Ве- 
личина дискрет быстро возрастает с увеличением к, и 
это увеличение тем значительнее, чем больше значение 
к. В связи с быстрым нарастанием величины дискрет 
дифференциального спектра возникает проблема пере- 
полнения разрядной сетки при компьютерном расчете 
дискрет. С одной стороны, для более точных вычис- 
лений требуется нахождение большого количества 
дискрет, а с другой стороны, при реальных величинах 
коэффициентов в дифференциальных уравнениях при 
вычислении дискрет с большим номером к неизбежно 
появляются ошибки вычислений. В силу рекуррентного 
характера вычисления дискрет ошибка вычисления 
накапливается. 

Цель работы - исследование возможности разреше- 
ния этого противоречия путем перехода к дифференци- 
ально-чебышевскому базису. 

Рассмотрим переход к базису ортогональных мно- 
гочленов Чебышева первого рода и смещенных много- 
членов Чебышева. Для этого необходимо установить 
связь между свободными коэффициентами степенного 
и анализируемого базисов. Результат может считаться 
положительным, если при переходе к полиномиальному 
базису скорость сходимости ряда увеличится. 


2. Переход к дифференциально-чебышевскому базису 

Среди всех полиномиальных базисов наилучшими 
свойствами сходимости обладают многочлены Чебыше- 
ва [2] первого рода ЦД), образующие базис {Т к (Ц} . Они 
обеспечивают наилучшую аппроксимацию в смысле 
минимаксной нормы [3]. Переход к чебышевскому 

леной 
мно- 
гочленами Чебышева. 

Формула перехода в базис ортогональных полино- 
мов Чебышева первого рода Т^І) получена в [2] и может 
быть записана в виде: 
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Т ; (1) = еок(1 -агссо.ч(1)) 7 (3) 

где т=0; п - номер старшего многочлена Чебышева, а 
{.} - целая часть числа. 

В качестве примера рассмотрим функцию 

? 0 О=ьа+ т )«Е(-і) ,+1т ’/ і = ? 1 ( т у 

м 

для которой с помощью формулы (3) получаем следу- 
ющее разложение: 

^ 2 О = -0,39б т 0 О + 1 ; 375 т 1 О- 

-0,453 Т 2 ( т ) + 0,146 Т 3 ( т )- < 4 ) 

-0,0б2 Т 4 О+0,013 Т 3 О-0,005 Т 6 О. 


В разложении (4) отбрасывание последнего слага- 
емого приводит к ошибке 8=0,17, а в разложении (5) 
отбрасывание трех последних слагаемых приведет к 
ошибке 8=0,09. Построим для сравнения две зависи- 
мости, которые представляют собой полиномы третьей 
степени. Первый полином (6) образован из исходной 
(4) отбрасыванием членов выше третьего порядка, а 
второй (7) представляет собой полином так же третьей 
степени, но полученный с помощью перехода к базису 
разложения по первого рода. 

х,(о= інгѵ/і, 

І=1 

^ 4 О = -0,396 Т 0 О + 1,375 Т 1 ( т )- 
-0,453 Т 2 ( т ) + 0,146 Т 3 ( т ), 


где Т 0 (1)=1, Т^І)^, Т 2 (1)=2 і 2 -1, Т 3 (1)=4і 3 -Зі - многочле- 
ны Чебышева первого рода. Окончательно получаем 

Ч О = 0, 057 + 0, 93 7 Т - 0, 906° + 0,584 т3 . (6) 

Из зависимостей, построенных по формулам (6) и 
(7), наглядно видно, что при одинаковом порядке по- 
линомов аппроксимация функции х(1) точнее в случае 
использования чебышевского базиса (см. рис. 1). В этом 
случае максимальное отклонение х 4 (1) от аппроксими- 
руемой функции х 2 (1) имеет меньшее значение. 

Вышеизложенную методику перехода к чебышев- 
скому базису можно применить, в частности, в т® 
случае, когда имеется некоторый полином, являющийся 
аппроксимацией неизвестной функции, и требуется по- 
низить порядок этого полинома. Такая ситуация может 
возникнуть при экспериментальных исследованиях, 
результатом которых являются значения, являющиеся 
отсчетами какого-либо процесса. 

Если вернуться к разложению (3), можно увидеть, 
что последние два коэффициента чебышевского базиса 
зависят только от одного соответствующего коэффици- 
ента степенного базиса. Таким образом, для вычисления 


двух последних дискрет чебышевского базиса формулу 
(3) можно преобразовать к следующему виду: 
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То есть, чем больше количество учитываемых диск- 
рет, тем значительнее будет уменьшение их величины в 
чебышевском базисе относительно степенного. Причем 
отношение двух последних дискрет при Н=1 всегда 
будет равно соответствующей степени числа два. Если 
же Н будет больше единицы, дискреты чебышевского 
спектра будут уменьшаться как соответствующая сте- 
пень удвоенного Н, то есть еще более быстро. Таким 
образом, показано, что в чебышевском базисе дискреты 
спектра убывают значительно быстрее, чем в степенном, 
а скорость сходимости ряда увеличивается. /о\ 


3. Дифференциально-чебышевские преобразования со смещен- 
ными полиномами Чебышева 

Авторы [2] упустили возможность исследования 
точности вычисления дискрет в базисе смещенных 
полиномов Чебышева [3]. Проведем эти преобразова- 
ния. Значение смещенных полиномов Чебышева 8^1) 
вычисляются по следующей формуле: 
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Первые шесть смещенных многочленов Чебышева 
имеют вид: 

$о(9 = 1; 

3,(1) = 21-1; 

8 2 (1) = 8і 2 -8і+1; 

8 3 (I) = 32т 3 - 48і 2 +181-1; 

§ 4 (I) = 128Т 4 - 256т 3 + 1 60( 2 - 321 + 1 ; 

8 5 (I) = 5 12(5 - 1280і 4 +1 120Т 3 - 400т 2 + 501-1 ; 

8 6 (I) = 2048Т 6 - 6144т 5 + 6912т 4 — 3584т 3 + 840т 2 - 721+1. 

Отсюда последовательно находятся степени отно- 
сительного аргумента ()=і/Н как линейные комбинату 
многочленов Чебышева: 
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® 3 =[Ю8 0 ( ® )+158,( ® )+68 2 ( ® Ч8 3 ( Ѳ г ;32. 

® 4 =[358 0 ( ® )+568,(® )+288 2 ( ® )+88 3 (® )+3 4 (® )]/128, 
ит.д. 

Обозначая коэффициенты данных зависимостей 
через С і? можно 

ѳ' =Х х * І *(Ѳ) 

1 =0 
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Отсюда оригинал х(1), имеющий дифференциаль- 
ный спектр Х(к), теперь может быть представлен в 
виде ряда: 
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Рис. 1. Аппроксимация функции 1п(1+і) 


Ѵ)=Ё х ло 

=0 

Таким образом, связь коэффициентов разложения 
в степенном базисе и коэффициентов разложения 
в базисе смещенных многочленов Чебышева {8 і (і)} 
представляет собой формулы (9-11). 

Преимущество, которое получается при использова- 
нии разложения по смещенным многочленам Чебышева 
в смысле скорости сходимости ряда перед разложением 
по многочленам Чебышева первого рода, связано с тем, 
что в записи смещенных многочленов присутствуют все 
степени 1, а в записи многочленов первого рода іфй) 
сутствуют только либо четные, либо нечетные степени. 
То есть в первом случае коррекции подвергаются все 
члены полученного степенного ряда, а во втором лишь 
половина. Покажем это на примере с той же функцией 
1п( 1 +1) при п=6. Произведя разложение по смещенным 
полиномам Чебышева, получим 

^ 5 О = 0,34 і 0 О + 0,314 і 1 О- 

-0,048 е 2 О-5,53-1СГ зЕ зО- 
-3, 42 ■ ІО' 3 2 4 О - 5,86 -1 0 -4 % О - 
-8,14-10 5І 6 ( Т ). 

При сравнении полученной зависимости с выраже- 
нием (7) видно, что величина коэффициентов разложе- 
ния по смещенным многочленам Чебышева убывает 
значительно быстрее, что говорит о более высокой 
скорости сходимости ряда. 

В формуле перехода к коэффициентам разложения в 
чебышевском базисе (3) значение Н присутствует, в то 
время как в формуле (10) значение масштабной посто- 
янной не учитывается. Получается, что коэффициенты 
разложения по многочленам Чебышева первого рода за- 
висят от Н, а коэффициенты разложения по смещенным 
многочленам Чебышева при изменении Н постоянны. 


Значение Н учитывается при подстановке выражений 
многочленов для получения решения в виде степенного 
ряда, то есть на самом последнем этапе. 

Если вычислить по формуле (10) последний коэф- 
фициент разложения в чебышевском базисе, то он будет 
зависеть только от одного коэффициента разложения в 
степенном базисе. Тогда для вычисления последнего 
коэффициента разложения по смещенным многочленам 
Чебышева формула (101 шэеобшзуется к виду 

г _ Х(п) 

п 2 2п - ] * 

При сравнении формулы (2) с формулой (8) для 
последнего коэффициента разложения по многочленам 
Чебышева первого рода видно, что величина последней 
дискреты в баз тт ^ л многочленов Чебышева первого рода 
будет всегда в ™ раз больше, чем величина последней 
дискреты в базисе смещенных многочленов Чебышева, 
так как 

Х(п) 2 м = 2 п 
Г 1 ' Х(п) 

Суть того, что лучшая аппроксимация функции 
получается при отбрасывании членов ряда не в степен- 
ном базисе, а в базисах ортогональных многочлен^) 
заключается в следующем. При отбрасывании члена 
ряда максимального порядка в степенном базисе осталь- 
ное выражение остается без изменения, в то время как 
при отбрасывании члена ряда максимального порядка 
в базисах ортогональных многочленов за счет их вида 
происходит коррекция коэффициентов степенного ряда, 
образующегося при подстановке соответствующих ба- 
зисных многочленов. То есть при понижении степени 
функции автоматически корректируются коэффициенты 
при меньших степенях. 

4. Решение плохообусловленных систем уравнений с помощью 

перехода к чебышевским базисам 
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и результат восстановления при 10, 13 и 20 значениях 
функции, выбранных с равномерным шагом. Видно, 
что при увеличении числа отсчетов аппроксимирующая 
функция приближается к истинной. 

Переход к разложению по многочленам Чебышева 
можно использовать не только при работе с дифферен- 
циальным преобразованием, но и при решении плохо- 
обусловленных систем больших порядков. Для этого 
необходимо решить систему уравнений в чебышевском 
базисе, а затем подставить соответствующие многочле- 
ны Чебышева в полученное выражение. 


Рис. 2. 


Аппроксимация функции 
с ростом числа отсчетов 


1(1) = 2,Зехр(4()5Іп(61) 


В радиофизике приходится часто сталкиваться с 
решением плохообусловленных систем уравнений. 
Например, при исследовании временных процессов 
часто требуется найти достаточно точное аппроксими- 
рующее выражение, что связано с решением как раз 
плохообусловленных систем уравнений, так как это 
довольно кратковременные процессы. Поэтому для 
решения плохообусловленной системы уравнений, не 
разрешимой в степенном базисе, необходим предло- 
женный переход к чебышевскому базису: 


^а і Т і Ц і )=ТЦ і ), і = 1,п; і=1,п, 

і=0 


где Т п (1 п ) - численное значение многочлена Чебышева 
при соответствующем значении 1; Т(1 П ) - значение фун- 
кции при соответствующем значении 1; а п - искомые 
коэффициенты. Записанная таким образом система 
лишена недостатков, характерных для плохо обуслов- 
ленных систем. 




тг о тт а л т п а тт п ті лт а 


ра возьмем функцию 
Г(Е) - 2, 3 ехр(41) ып(Ѳ1:). у схан0 влено, что в сте- 
пенном базисе максимальное количество отсчетов, 
по которому можно восстановить функцию по ур. (2), 
равно десяти для таких вычислительных пакетов, как 
МаИіСАБ и МаіЫаЪ. Для большего количества зна- 
чений ошибка восстановления слишком велика, или 
итерационный процесс вообще расходится. С помощью 
перехода к разложению по многочленам Чебышева пер- 
вого рода можно найти решение при гораздо большем 
числе отсчетов. На рис. 2 изображена исходная функция 


5. Заключение 

Таким образом, в работе показано, что с помощью 
перехода к чебышевским полиномиальным базисам 
увеличивается скорость сходимости ряда, в виде ко- 
торого получается решение уравнения. Разложение по 
смещенным многочленам Чебышева наиболее предпоч- 
тительно. В чебышевских базисах можно оценивать точ- 
ность полученного решения, последовательно вычисляя 
дискреты, так как величина чебышевских дискрет 
постоянно уменьшается с увеличением их номера. Так 
же показано, что с помощью перехода к чебышевским 
полиномиальным базисам появляется возможность 
решения систем больших порядков, решение которых 
обычным способом затруднительно. 
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